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ABSTRACT 



Of 



Recent evidence for a strong 'hard excess' of flux at energies > 20 keV in some Suzaku 
observations of type 1 Active Galactic Nuclei (AGN) has motivated an exploratory study 
of the phenomenon in the local type 1 AGN population. We have selected all type 1 AGN 
in the Swift Burst Alert Telescope (BAT) 58-month catalog and cross-correlated them 
with the holdings of the Suzaku public archive. We find the hard excess phenomenon 
to be a ubiquitous property of type 1 AGN. Taken together, the spectral hardness and 
equivalent width of Fe Ka emission are consistent with reprocessing by an ensemble of 
Compton-thick clouds that partially cover the continuum source. In the context of such 
a model, ~ 80% of the sample has a hardness ratio consistent with > 50% covering of 
the continuum by low-ionization, Compton-thick gas. More detailed study of the three 
hardest X-ray spectra in our sample reveal a sharp Fe K absorption edge at ~ 7 keV 
in each of them, indicating that blurred reflection is not responsible for the very hard 
spectral forms. Simple considerations place the distribution of Compton-thick clouds 
at or within the optical broad line region. 



Subject headings: galaxies: active - X-rays: galaxies - Seyfert - X-rays 



1. Introduction 

The spectral shape of the cosmic X-ray background suggests that a fraction of active galactic 
nuclei (AGN) are heavily obscured by Compton-thick mat erial; however, the significance of tha t 
contribution to the cosmic X-ray background is debatable (jGilli et al.1 120071 : iTreister et al.1 120091 ) . 
Nonetheless, the extensive search for Compton-thick AGN has lead some to identify these sources 
via the r atio between the X-ray flux and either the mid-infrared (MIR ) flux, [Oiii] line or [Ne v] 



line (e.g. iHeckman et al. 



2005 



Gilli et al 



2010; 



LaMassa et al 



2011 



quantities assumed to be 
negligibly affected by obscuration, and therefore representative of the AGN bolometric flux. Other 
groups have identified Compton-thick sources via the hardness ratios between the medium and soft 
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X-ray bands (e.g. IWinter et al.l 12009 ; iBrightman fc Nandrall2012l ). believing this bandpass to be 
best for the detection of X-ray absorption. 

The cosmic X-ray background is known to require AGN to have a distributi on of obscuring 



column densities that extend into the Com pton thick regime, Nh > 10 cm (c.f. IComastri et al 



19951 : iGilli et al.1 120071 : iTreister et al,ll2009l ). and it has been assumed that the most obscured AGN 
have hitherto largely escaped direct detection. However, there is now widespread evidence for 
partial covering of the X-ray absorber even i n type I AGN that in those models have previously 



been assum ed to be largely unobscured (e.g. iTanaka et al.l I2004J : iMiller et al.l 120071 : iTurner et al 
20091 . boilh . 



High resolution UV s pectroscopy has revealed mu l ti-layered, complex ab sorption to be a com- 
mon phenomenon in AGN (| Crenshaw &; Kraemerll 19991 : ICrenshaw et al.ll2003l ) . The det ection of ab- 



sorpt ion features, such as H- and He-like species of C, N, O, Ne, Mg, Al, Si, and S (e.g.. iKaspi et al 
2002), has shown that the signature s of complex absorp t ion extend into the X-r ay band, with warm 
absorbers being common in AGN (jBlustin et al.ll2005l : iMcKernan et al.ll2007l ). Initially, complex 



absorption in the soft X-ray band was detected over three decades in ionization p arameter (log 



£ ~ - 3) and three deca des in column density (N# 



10 



20 



1Q 23 



cm 



2003 



e.g. 



Netzer et al 



Turner fe Milled 120091 ). However, the detection of deep Fe XXV and Fe XXVI absorption lines 



expanded the known range of ionization parameter and column density in t he AGN population . 



Observati ons of NGC 3783 (jKaspi et al 



2002 



Mrk 766 (IMiller et al.ll2007ft and NGC 3516 (ITurner et al 



Reeves et a 



2004), NGC 1365 (jBisaliti et al.ll2005h 



20081 ) were among the first to show deep 



K-shell absorption lines from highly-ionized species of Fe. 



Recently, iTombesi et al.l (|2010i ) performed a blind search for absorption signatures in the Fe 
K regime of radio-quiet AGN by cross-correlating the sources detected in the RXTE All-Sky Slew 
Survey with the XMM-Newton archive and detected absorption from Fe xxv, Fe xxvi, or both, in 
~ 40% of the sources studied in their sample. The depth of these absorption lines constrains this 



highly ionized component of gas to iV#~ 10 - 10 cm along the line-of-sight. Thus, the full 
set of absorption signatures observed throughout the X-ray band have shown the X-ray absorbing 
gas to extend over a range of column density (N# ~ 10 20 - few x 10 24 cm" 2 ), ionization parameter 
(log £ ~ - 5), an d outflow velocity a few hundred km s _1 to 0.3 c (e.g. IMcKernan et al. 



2007 



Tombesi et al. 2011 



In addition to there being a wide range of gas conditions across the AGN p opulation, it 



is also evident that indiv id ual source s requ ire multiple zones of absorbing gas (e.g. iNetzer et al 



2003 



Blustin et al.1 120071 ). iLee et al.l (|200ll ) showed that MCG-6-30-15 required the presence of 



at least two distin ct zones of g as to explain the soft-band absorption features in the Chandra 
MEG grating data. IMiller et al.l (|2008l ) demonstrated that one could extend the complex absorber 
model to include zones of a higher column density, whose variations in covering fraction could 
then explain the marked spectral variability observed in MCG-6-30-15, with no requirement for 
a contribution from blurred reflection. This multi-epoch analysis revealed MCG-6-30-15 to be 
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typical for a general class of AGN dominated by absorption effects. Similarly, a study of NGC 
3516 found that the emission and absorption features detected in the grating data requi red the 
presence of four distinguishable zones of gas to shape the spectrum ([Turner et al.ll2008l ). The 
marked spectral variability in this source was attributed to variations in the covering fractions of 
one of the inter mediate layers. Notably, a broad dip has been observed in the X-ray light curve 



of both objects (jMcKernan &: Yaqoob 



1998 



Turner et al.ll2008l ). taken as supporting evidence for 



occultation of the X-ray continuum by absorbing gas. 

Further evidence for large columns of absorbing gas along the line-of-sight has come from 
Suzaku. The Suzaku observations of 1H 0419-577, PDS 456 and NGC 1365 revealed a marked 
excess of flux a bove 20 keV, compared to that predic t ed from fits to data below 10 keV, dubbed 



a 'hard excess' (Turner et al 



2009; 



Reeves et al 



2009 



Risaliti et al.1 l2009al . respectively) . In these 



sources, the high PIN-band flux (15-50 keV) was explained by the presence of a Compton-thick 
absorber covering > 70% of the continuum source, suggesting that type 1 AGN can be unobscured 
in the optical bandpass, but partially covered in the X-ray regime. In order to obtain the true 
intrinsic X-ray luminosity of these sou rces, one must app l y significant correc tions for absorption 



and for Compton-scattering losses (e.g. iTurner et al. 



2009 



Reeves et al J 12003 ). 



The results from 1H 0419-577, PDS 456, and NGC 1365 motivated an exploratory study of 
the hard excess phenomenon in the local type 1 AGN population. In this paper, we investigate 
this phenomenon using a BAT-selected sample of type 1 AGN. We begin in Section 2 by outlining 
the observations and data reduction. In section 3, we measure spectral hardness ratios and FeKa 
emission line equivalent widths for the sample, and present a simple Monte Carlo radiative transfer 
calculation that can reproduce these measurements.. In Section 4, we present the results of some 
case studies. In Section 5, we discusses the implications of our analysis. We draw conclusions in 
Section 6. 



2. The Observational Data 

The Swift BAT is a sensitive coded aperture imaging instrument with a large (1.4 steradian) 
field-of-view. The energy range over which BAT can perform imaging observations is 15-150 keV, 
yielding positions accurate to ~ 4'. The bandpass of BAT allows it to make an unbiased survey 
of the X-ray sky for column densities up to ~ 10 24 cm -2 . We selected all type 1 AGN, including 
intermediates up to type 1.9 and exclud ing known radio sources and LINERS, from the 58- month 
BAT catalog ( Baumgartner et al. 2010 ). To understand the sample properties, we required simul- 



taneous medium (2-10 keV) and hard X-ray (>10 keV) data. Therefore, we cross-correlated those 
selected sources with the Suzaku archive, including all sources available in the public domain as of 
2011 December 20. Our selection criteria yielded a total sample of 43 objects and 76 observations: 



http: / /heasarc. gsfc.nasa.gov / docs/swift / results/bs58mon / 
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24 objects are classified as type 1-1.2, 16 objects are classified as type 1.5, and 3 objects are clas- 
sified as type 1.8-1.9. Sources with multiple observations were reduced and analyzed individually. 
Details concerning the observations and data are presented in Table [TJ 

Suzaku has four X-ray telescopes, each containing a silicon CCD within its focal plane forming 
the X-ray Imaging Spectrometers (XIS) suite. XISO, XIS2 and XIS3 are front- illuminated (FI), 
providing data over a usable range of 0.6-10.0 keV with an energy resolution FWHM ~ 130 eV at 
6.0 keV. In November 2006, a charge leak was discovered in XIS2, making XISO and XIS3 the only 
operational FI chips. XIS1 is back-illuminated. The back- illuminated configuration extends the 
soft band to ~0.2 keV; however, it also results in a lower effective area and higher background rate 
in the Fe K regime, compared to the FI chips. Consequently, XIS1 is excluded from our spectral 
analysis. Suzaku also carries the Hard X-ray Detector (HXD) that contains a silicon PIN diode 
detector covering a range of 10-100 keV with a usable energy range of 15-50 keV. 

The data were reduced using HEAsoft v. 6. 10. The XIS cleaned event files were screened in 
XSELECT to exclude data during passage through the South Atlantic Anomaly and also excluding 
data starting 500 s before entry and up to 500 s after exit. In addition, we excluded data having 
an Earth elevation angle < 10° and a cut-off rigidity > 6 GeV. CCDs were in 3 x 3 and 5x5 edit 
modes, with normal clocking mode. Good events were selected, having grades 0, 2, 3, 4, and 6, 
while hot and flickering pixels were removed using the SIS CLEAN script. The spaced-row charge 
injection was utilized. XIS spectra were extracted from circular regions of 3.0' radius centered on 
the source, while the background was extracted from a region of the same size offset from the source 
and from the corners of the chip that register calibration data. 

The cleaned PIN data was reduced utilizing the ftool hxdpinxbpi. This tool calculates good 
time intervals (GTIs) of the non X-ray instrumental background (NXB, using model 'D' released 
2008 June 170 ) data that overlap with the source data and extracts both the source spectrum 
and NXB background spectrum through that common GTI. A simulated cosmic X-ray background 



(CXB) spectrum is calculated, based on the CXB spectrum reported in iBoldtl (|1987l ) and renor- 
malized to the 35' x 35' field of view of the PIN instrument, and combined with the NXB spectrum 
to produce a total PIN background spectrum for the observation. A dead-time correction ~ 4-5% 
is applied to the source spectrum using the ftool hxddtcor. Finally, the PIN data have a known 
la systematic uncertainty of 1.39cu, which was applied to the PIN data in grppha. 



3. Spectral Analysis and Sample Results 

In this paper we are primarily interested in the gross spectral properties of the sample of 
AGN, and in particular in the broad-band hardness ratios and emission line equivalent widths. In 



2 http://www.astro. isas.jaxa.jp/suzaku/doc/suzakumemo/suzakumemo-2007-01.pdf 
3 http: / /heasarc. nasa. gov / docs /suzaku / analysis / watchout . html 
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order to measure these quantities, we fit some simple models to each AGN spectrum. The aim of 
these models is not to obtain accurate physical parameters for the X-ray sources, but rather to get 
sufficiently accurate measurements of broad-band fluxes, corrected for the instrumental response, 
and sufficiently accurate measurements of FeKa emission line equivalent width, to be able to 
characterize the gross properties of the spectra. 

We fitted each observation over the 2-50 ke V bandpass using XSPEC v 12.5 . The model 
used was a powerlaw modified by an xstar grid (jKallman Bautistall200ll ) that was allowed to 
partially-cover the continuum. Further to this, we allowed a Gaussian component at ~ 6.4 keV. 
Finally, a column of gas representing the Galactic line-of-sight absorption was parameterized using 
the tbabs model ( Wilms et alj|2000h and w as fixed to the weighted average N# in the Dickey and 
Lockman survey ([Dickey fe Lockmanlll990l ). The XSTAR table was based on v2.2.0 of the code and 
had a gas density n = 10 10 cm -3 , an illuminating photon index F = 2.2 and a gas turbulent velocity 
300 km s _1 . In AGN, the gas density is poorly constrained and is degenerate with the ionization 
parameter and assumed radial location of the gas. As xstar requires the gas density as an input 
parameter, we chose a value that lies at a point where collision effects and continuum dilution are 
negligible. The photon index is justified in Section 3.3.2; the turbulent velocity was based on the 
approximate outflow velocities of the warm absorbers. 

During spectral fitting the normalization and photon index, T, of the powerlaw continuum and 
the column density, covering fraction and ionization parameter of the xstar table were allowed to 
vary. This model provided an adequate representation of the observed continuum, from which we 
could determine the flux in each band of interest in a model independent way. The mean reduced 
X 2 value for the sample when fitted to our model was ~1.5. Ten observations had reduced \ 2 > 
2. More complex models yielded a reduced x 2 ~ 1 f° r these 10 observations, and the difference in 
the hardness ratio calculated from the simple model and the more complex model was < 15% for 
each observation. We concluded that our simple model fits gave an accurate measurement of the 
integrated flux in the bands of interest. 

During spectral fitting, the PIN part of the model was scaled by an energy-independent con- 
stant as appropriate for the XIS or HXD nominal aim point used, 1.16 or 1.18, respectively. This 
scaling factor accounts for the cross-calibration constant required to correct ly reconcile the XIS and 
PIN data based on the calibration of the Crab Nebula (jMaeda et al.l l2008). 



We extracted the observed energy fluxes (ergs cm -2 s~ x ) for the 2-10 keV and 15-50 keV 
bandpasses to determine the hardness ratio, Fluxi5_5o keV /Flux2_io keV, fo r each observation. These 
bandpasses were chosen as the most meaningful and practical for determination of the hardness of 
the X-ray spectrum associated with the active nucleus. In this field, the 2-10 keV band has become a 
standard bandpass over which to quote a flux, this is partly because many previous X-ray detectors 
covered this bandpass (e.g. the EXOSAT ME, ASCA SIS and GIS, X MM EPIC CCDs). In any 
case, below 2 keV there are significant contributions to the X-ray flux from extended gas, starburst 
regions and other structures that are not immediately associated with the nucleus. There are no 
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reliable data available between 10-15 keV. On the high end, the upper limit was chosen based on 
the limited sensitivity of the PIN above 50 keV. 

To obtain the intrinsic hardness of the AGN, one should correct for the Galactic absorption. 
However, we found that the Galactic absorption had a negligible effect on the 2-10 keV and 15-50 
keV flux measurements for the entire sample. 



3.1. The X-ray Hardness 



Figure Q] shows the hardness rati o, Flux^^n^y/Flu x? 



-inkeVi f° r each observation in our 
sample, plotted against the BAT flux (jBaumgartner et al.ll2010i ). The black line is the weighted 



mean hardness ratio (1.65 ± 0.01) of the sample. 

To understand the distribution of hardness ratios in the sample, we compared our results to 
some simple models. For these model s, we assumed r= 2.0, consistent with the average photon index 
found for a sample of type 1 AGN (cf. lScott et al.ll201ll ) and a continuum cut-off at 500 keV, beyond 
the bandpass considered here. First we calculated the expected ha rdness ratio for a simple, dis k 
reflection model. The reflection was parameterized using pexrav (jMagdziarz fc Zdziarskilll995l ). 
and represents that from a standard thin disk of neutral material subtending 2ir sr at the continuum 
source and having Solar abundances. The disk was assumed to be observed at an inclination of 
60°, as the model is insensitive to inclination (showing only a ~ 10% difference in the hardness 
ratio relative to our chosen value for inclinations of 30° and 80°. The model hardness ratio shown 
is for the sum of the powerlaw and reflection. In our sample, 90% of the objects are harder than 
this model prediction. 

We also compared our observational result to the model prediction where the continuum source 
is completely hidden (pure reflection) , using the same parameter values as the previous disk model. 
Interestingly, even pure reflection cannot explain the most extreme objects in this sample (Fig. 1). 
We note that the sources above the pure reflection model line are typically type 1.8 and 1.9. While 
pure reflection may be expected from type 2 AGN, such an explanation is at odds with a sample 
of type 1 AGN in the context of the current paradigm. 

Of course, the hardness ratio for this model (and the others discussed below) depends on the 
photon inde x assumed. From theory, one might expect the photon index to lie in the range 1.8< 
r < 2.5 (cf. lHaardt &: Maraschilll993l ): for the simple reflection model noted above, the hardness 
ratio would be 1.5 for an illuminating continuum of photon index T = 1.8, 1.1 for T = 2.0, and 0.4 
for r = 2.5. 

We also compared the data to predictions based upon a simple partial-covering model, char- 
acterized by a neutral absorber (pcfabs in xspec). This simple model does not include Compton- 
scattering as such an effect is geometry dependent. Including Compton-scattering requires an 
understanding of the distribution and location of the Compton-thick gas, which currently is not 
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well constrained for type 1 AGN. Markers on the right side of Figure [T] represent the ratios for the 
case where a neutral column of Compton-thick gas partially covers the T = 2 continuum source, 
taken here to be 2 x 10 24 cm -2 . Several covering fractions are denoted: 98% (green), 90% (blue), 
70% (cyan) and 50% (magenta). Approximately 80% (35 objects, 63 observations) of the sample 
have hardness ratios consistent with > 50% covering of the continuum by Compton-thick gas. 



3.2. The Equivalent Width of Fe Ka Emission 



The broad-band X-ray spectra of local Seyfert type 1 AGN are much harder than was expected, 
based on an excess in the observed flux above 20 keV, compared to that predicted from fits to 
data below 10 keV. From a simple consideration of the distribution of hardness ratios we have 
demonstrated that these X-ray spectra are shaped by a high covering fraction of neutral gas along 
the line-of-sight, or are dominated by a reflected X-ray component. 

Inspection of the Fe Ka line properties can clearly help us to distinguish the origin of spectral 
hardness. For example, if we were observing the reflected component without seeing the illuminating 
continuum, then we would e xpect a large (> 1 keV) l ine equivalent width (EW), measured agains t 



the reflected continuum (e.g. iKrolik fe Kallman 



1987 



George fe Fabian 



1991 



LaMassa et al 



2011 



In contrast, partial covering absorber models predict smaller line equivalent widths when measured 
against the observed continuum (e.g. iMiller et alj|2009l : lYaqoob et al.ll2010i ). 



Suzaku has a high throughput in the Fe K regime, allowing for accurate constraints to be placed 
on the Fe Ka emission line for our sample sources. The total Fe Ka emission line was parameterized 
using a Gaussian profile with a freely varying width (<r), normalization and energy. In 9 sources 
(11 observations), the statistical fit improved significantly (A% 2 > 10 and F-test probability < 0.05 
) when an additional Gaussian line component, to parameterize a broader component of Fe Ka 
emission, was allowed in the fit (with all parameters allowed to be free). In these sources, the two 
Gaussian components accounted for the core and broad base in the spectral shape of the total Fe Ka 
emission. We determined the equivalent width of the total Fe Ka emission by calculating it against 
the total observed continuum. The average total equivalent width for those observations with one 
Gaussian component and for those observations with narrow and broad Gaussian components were 
both ~ 150 eV. 

Additionally, 10 observations showed significant (Ax 2 > 10) line emission from Fe xxv, Fe 
XXVI, or both. These features were parameterized with a narrow (a = 10 eV) Gaussian component 
fixed at 6.7 keV or 6.97 keV, respectively. We will present the results from our detailed spectra l 
analysis, including the EW of emission in the Fe K regime, in our next paper (jTatum et al.ll2012al ). 

Figure [2ji shows the hardness ratio plotted against the EW of the total Fe Ka emission line. 
Several models are overlaid on the data, and for the models considered, the Fe Ka emission line 
EW has a dependence on photon index, i nclination and iron abundance of the reprocessor (c.f. 



George fe Fabian 



1991 



Nandra et al.l 120071 ). As before, we assumed V = 2.0 and solar abundances 
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for all models. 

First, we considered the spectral hardness and line equivalent width that might be produced 
from simple reflection of a power-law continuum from a standard accretion disk (as described 
previously). Figure 2a shows model predictions, parameterized using pexrav, for R=l (against the 
illuminating continuum) through to pure reflection (i? — > oo). It is clear that for a given hardness 
ratio, disk reflection models consistently predict a much larger Fe line EW than is observed. 



We also compared our EW distribution to the values predicted by MYTORUS (jMurphy k, Yaqoob 



2009), a torodial reprocessor valid in the Compton-thick regime. In our MYTorus tables, angles 
in the range 0° < 9 < 60° represent an approximately face-on view, where the continuum is not 
obscured by the toroidal material. As the MYTorus spectra have only a very weak dependance on 
inclination angle within this range (with regard to our particular quantities for comparison), for 
clarity, we show only the model line for 6 = 60° as a representation of the 'face-on' view. The 
small loop (Fig. 2a) shows the MYTorus predictions as the column density of the torodial gas is 
stepped over the range 10 23 - 10 25 cm . For a face-on torus, the hardness ratio peaks at ~ 1.0, 
falling below the bulk of the data points, while line equivalent widths (18 - 65 eV) do not achieve 
the span evident in the observed values. Exploring a more edge-on orientation for MYTorus, where 
the continuum is obscured, we fixed the inclination angle to 70° and varied the column density of 
the torodial gas as before. As illustrated in Figure 2a, such a case does produce model predictions 
that span the range of values observed. Interestingly, the range of inclinations between the face-on 
scenario and 70° would be consistent with the data points. However, if a large fraction of type 1 
AGN require large inclination angles for the obscuring torus then this presents a challenge for the 
Unification Model; the Seyfert 1 galaxies should represent the face-on cases. 

Whi le the possibility o f toroidal obscuration is intriguing, well-studied sources, such as MCG- 



6-30-15 ([Miller et al.ll2008l ). appear to have a spectral form shaped by a complex absorber, with 
some zones only partially covering the continuum source. X-ray spectral variability can then be 
explained as changes in the covering fraction of the absorbing gas. Such a model obviously requires 
the gas to be clumpy, rather than a solid and uniform form and indeed, it may be the case that 
there is simply a toroidal distribution of clouds in these AGN, producing observed properties close 
to the values for the MYTorus model. Another possibility is that changes in the continuum source 
size could appear to mimic covering changes, as the fractions of obscured and unobscured light 
would vary. However, such a model would place tight constraints on the continuum size and torus 
location. 



3.3. A Cloud Model for the Absorbing Gas 

The hardness ratio, Fe Ka line equivalent widths, the detailed X-ray spectral shape and vari- 
ability all point to a distribution of clouds, partially covering the X-ray source. To investigate this 
in more detail, we performed a Monte-carlo calculation of the X-ray spectrum expected from a dis- 
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tribution of cold gas clou ds surrounding a central source of X-rays. This idea was first explored by 



Nandra fc George! (| 19941 ) but computing power available at that time limited the ability to follow 



virtual photons as they scatter between clouds. Now, it is possible to follow photons as they scatter 
any number of times between multiple clouds. 



3.3.1. The distribution of scattering clouds 

In the calculation, gas is placed in a distribution surrounding a central X-ray source. The 
distribution is assumed to be a strict two-phase medium which is either vacuum or gas. The gas 
phase has a constant density everywhere. The distribution is created in two steps: 

1. A number of points are randomly placed with a uniform spatial density within a spherical 
shell centered on the primary source, with a defined inner and outer radius. 

2. The gas distribution is assumed to consist of spheres of constant gas density drawn around 
that set of points. If two (or more) spheres overlap, the density is not doubled, but remains 
at the same value as elsewhere in the gas phase. Thus in the limit where the density of sphere 
centers is very high, the distribution tends towards being a constant density spherical shell. 

For distributions with a high density of points, the spheres heavily overlap, and the resulting gas 
distribution does not then resemble a set of individual blobs, but has a sponge-like topology. For 
the radiative transfer calculation, the spherical distribution is assumed to be bisected by an infinite, 
thin accretion disk, that absorbs photons but does not itself radiate at X-ray energies, so that the 
scattered light received by the observer is reduced by a factor two compared with the optically-thin 
case without the disk. We dub the calculation "Monte-Carlo Radiative Transfer model (MCRT)". 



3.3.2. The opacity, line radiation and scattering cross- section 



The primary source of photons is assumed to be a point-like source producing a power-law X- 
ray spectrum with photon index V = 2.2, consistent with the photon index found fo r the best studied 
sources once all the absorbi ng zones wer e taken into accoun t (i.e., MCG-6-30-15 iMiller et alJl200a . 
Mrk 766 lTurner et al.ll2007l . lH0419-577 lTurner et aHhood ). The gas is assumed to be cool, so that 
all ions heavier than H are neutral, but with H ionized, so that the number of free scattering electrons 
equals the number of H ato ms. The absorption o pacity is considered to be due to all elements with 
Solar abundances given by lAnders fc Grevessd (|1989l ). In reality the Compton-scattering clouds 
are likely to be hot and ionized, however a full treatment of ra diative transfer through such gas 
is computationally extremely expensive (see e.g. ISim et al.ll2010l ) and for the purpose of exploring 
a wide range of parameter space, the model presente d here makes the simplifying assumption of 
cold gas clouds. In a future work (jTatum et al.ll2012al ). we will compare our results with those that 
consider a full treatment of radiative transfer through a hot and ionized gas. 
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As photons propagate through the distribution they may be either absorbed by the gas, or 
be Compton scattered. After absorption, there may be resulting fluorescent line emission. This is 
implemented by calculating the mean free path for the photon in the gas due to three processes: 

1. Compton scattering, using the full energy-dependent Klein-Nishina cross-section, integrated 
over all scattering angles. 

2. Absorption by elements other than K-shell transitions of Fe, in which case the photon disap- 
pears from the calculation (i.e. line emission from elements other than Fe K-shell transitions 
is neglected). 

3. Absorption by an Fe K-shell transition, which may be followed by FeK fluorescent line emis- 
sion. 

3.3.3. Results from MCRT 

MCRT is a representation of neutral gas clouds partially covering an X-ray continuum source. 
Our calculations were performed for the case of 1000 cloud centers distributed within the half a 
spherical shell visible to the observer. 

The inner radius of the distribution is fixed at 10 units, scaled to the radius of an individual 
sphere. The outer radii have values 12.5, 14, 15, 17.5, 20 units, from the outer loop through to the 
innermost loop, respectively, corresponding to leaking light fractions of 50% - 70% in the 15-50 keV 
band. Distributions with higher volume filling factors have higher hard excesses and line equivalent 
widths. We calculated the hardness ratio/EW for a range of column densities covering 

10 23_ 10 25 

cm -2 for each of the filling factors considered. The model loops connect the model predicted values 
for different column densities, for a given filling factor. Moving around the loop corresponds to an 
individual sphere having a larger column density: on the outer curve are marked sphere column 
densities of logio(N^/cm _2 )= 24, 24.5, 25.0. As column density increases, photoelectric absorption 
hardens the spectrum and increases the fluorescent line strength. At high column densities, however, 
both the hard X-rays and the line flux from the high-column regions become highly attenuated, and 
the spectrum becomes dominated by a combination of the unobscured and scattered light, with an 
overall spectrum closer to the intrinsic spectrum of the source. For clarity, we plotted the model 
line predictions for this scenario separated from the data (Figure [2]). 

The looping model lines (solid red, dash green, dash-dot blue, dot cyan, dash-dot-dot-dot 
magenta) in Figure [2] are model predictions from MCRT with an increasing volume filling factor 
(traveling from the red curve to the magenta curve). The hardness ratios and Fe Ka emission line 
widths are well characterized by varying the volume filling factor of the neutral, clumpy gas, as the 
model lines lie within the la errors of 83% of the data. When only accounting for the narrow core 
of Fe Ka emission, the model lines would lie within the la errors of 85% of the data, suggesting 
that the narrow core comprises most of the Fe Ka emission. 
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MCRT can also calculate the expected relation for a thin, uniform absorbing shell, neglecting 
Compton scattering. This is shown in Figure [21 denoted as an orange model line. Moving from 
left to right along the curve corresponds to an increase in column density of the individual spheres 
ranging from 10 23 -10 25 cm . This model does not encompass the range of the observational data. 

Note that, in the absence of any other effects such as Doppler shifts, the results are independent 
of any scale size for the cloud distribution: all length scales are relative to the sphere radius, and the 
only dimensional parameter that enters the problem is the column density through an individual 
spherical cloud. 

Consideration of different filling factors of the clumpy gas show the MCRT to be consistent 
with the spectral hardness and line equivalent widths of the sample sources Figure El 



4. Case Studies 



We performed more extensive spectral analysis on the hardest sources and on a typical source 
to investigate the nature of the Compton-thick gas in more detail. 



4.1. A typical source: 1H 0419-577 

1H 0419-577 yielded one of the first detections of a hard excess i n a type 1 AGN. Th is source 



has been observed twice by Suzaku, (2007 July 25, 2010 January 16). iTurner et al.1 (J2009J) found a 
satisfactory fit to the 2007 July 25 observation of 1H 0419-577 using a partial covering model with 
16% and 67% of the line-of-sight obscured with N# ~ 5 x 10 23 and ~ 2 x 10 24 cm -2 and log £ ~ - 
0.1 and ~ 2, respectively (covering fractions not corrected for Compton-scattering losses). Analysis 
of the 2010 January 16 data revealed that the spectral shape remained unchanged between the two 
epochs, with the flux state of the source during the 2010 observation being 78% of that observed 
during 2007. This result verifies the existence of the detected hard excess and establishes it as a 
persistent feature of 1H 0419-577. 

Examination of the sample sources showed the surprising result that the hardness ratio (1.53 ± 
0.07) measured for 1H 0419-577 was consistent with the weighted mean of the sample (1.65 ± 0.01). 
Therefore, in the context of partial covering, the solution found for 1H 0419-577 could represent 
typical absorber parameters that may explain the hard excess phenomenon. 



4.2. The hardest objects 



The most extreme objects may provide the tightest constraint on the models. One naturally 
expects contributions from partial-covering absorption and accompanying reflection, in the context 
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of the MCRT, and together, these give a very hard spectral form, with spectral variability expected 
as the absorbed fraction changes. However, the alternative suggestion that might be invoked to 
explain our observational result is that the hard excess is the Compton hump of relativistically 
blurred disk reflecti on, with the core of the Fe line being reprocessed light from more distant gas. 



Walton et al.1 (|2010T ) applied the blurred reflection model to 1H 0419 - 577, PDS 456 and NGC 1365 



(2008 January 1). 

NGC 1194 (Sey 1.9) , NGC 1365 (Sey 1.8; X-ray data from 2010 July 15) and MCG 3-34- 
64 are the three observational datasets in our sample showing the largest hardness ratios. These 
datasets show a sharp Fe K absorption edge at ~ 7 keV (Figure [3]). To determine whether the 
edges were significantly blurred, we fitted the 5-9 keV band of each with a powerlaw modified by the 
Galactic column. We parameterized the 6.4 keV Fe Ka emission line with a Gaussian component 
and m odeled the Fe K abs orption edge near 7 keV with a smeared edge component (smedge in 
xsPEC. lEbisawa et al.lll994l ). All the parameters were allowed to vary, including the edge threshold 
energy, although this was not allowed to fall below the threshold energy of 7.11 keV for neutral Fe 
in the rest-frame. At 90% confidence, the depth and width of the edge for NGC 1194, NGC 1365 
(2010 July 15) and MCG 3-34-64 were r = 1.3 ± 0.40 and u < 0.15 keV, r = 0.87 ± 0.10 and 
a < 0.14 keV, and r = 0.95 ± 0.30 and a < 0.03 keV, respectively, showing that the absorption 
features are narrow and not likely to be reprocessed light from a blurred reflector. 

We performed an additional test for NGC 1365 as it had the highest quality data. We fitted 
the 6-9 keV bandpass using a powerlaw continuum with a Gaussian component for the Fe Ka 
emission, having a fixed intrinsic width a = lOeV. We then allowed blurring into the model, using 
the xspec model kdblur, with emissivity index q = 3, inclination angle 8 =,60° and the outer 
radius of the reflector at r ou t = 400r 9 . The data also required inclusion of absorption, which, in 
this limited bandpass, was adequately described by full covering from neutral gas with N# ~ 10 23 
cm . With this model, we found the best fit for the inner radius gave T{ n ~ 350 Tg. Stepping the 
inner radius down to rj n =20 r g found no other minima as the inner radius was decreased and the 
fit at rj n =20 r g was statistically worse by Ax 2 of 190 (F-test probability of ~ 10 -27 ), compared 
to the best fit. This test confirms that the reprocessor does not exhibit a blurred signature in this 
source. 

The most extreme sources possess sharp Fe K absorption edges. The absorption or reflection 
components that account for those edges also explain the hard spectral form of those sources. 
No relativistically-blurred components are required to explain these extreme sources, or any of the 
sample properties. We conclude that there is no evidence for any significant contribution by blurred 
reflection to the hard spectral form of local Seyfert galaxies. 
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Discussion 



The Suzaku observations of 1H 0419-577 (|Turner et al J 12003 ). PDS 456 (jReeves et all 12009 ) 



and NGC 1365 (jRisaliti et al.ll2009al ) were first to reveal a marked 'hard excess' of flux above 10 
keV. These results motivated an exploratory study of the hard excess phenomenon in the local type 
1 AGN population. 

We have selected all type 1 AGN in the BAT 58-month catalog and cross-correlated them 
with the holdings of the Suzaku public archive. We have found the hard excess phenomenon to 
be a ubiquitous property of type 1 AGN. Comparisons between the hardness ratio distribution 
and simple models suggested both neutral, partial-covering models and disk reflection models to 
be viable. Joint consideration of the spectral hardness and equivalent width of Fe Ka emission 
allowed us to rule out simple disk reflection models. 

The My Torus toroidal gas distribution accounts for the span of observed hardness ratios and 
line equivalent widths. However, evidence from detailed studies of nearby sources indicates that the 
reprocessor most likely is clumpy. We therefore performed simulations of an ensemble of Compton- 
thick reprocessing clouds. Predictions from the MCRT were in striking agreement with both the 
spectral hardness and the line equivalent width distributions. 

In the context of a simple absorption model assuming partial-covering by neutral gas, ~ 80% of 
the sample has hardness ratios consistent with > 50% covering of the continuum by Compton-thick 
material. Furthermore, the model predictions for the MCRT are within the la errors of 85% of 
the line equivalent widths in our sample. We conclude that both the spectral hardness and line 
equivalent width distribution strongly favor clumpy model interpretations. Our MCRT geometry 
shows excellent agreement with our observational results. 



5.1. Gas location 

5.1.1. Line width arguments 

Constraints obtained for the radial location of the absorbing gas are model dependent. Detailed 
analysis of the core of the Fe Ka emission line and spectral variability may place the t i ghtest 



constraints on the location of the circumnuclear gas in the context of the MCRT. iShu et al.l (|2010l . 
hereafter S10) studied the cores of the Fe Ka lines of 36 type 1 AGN using Chandra HEG (high- 
energy grating) data. The spectral resolution of HEG is unparalleled in the Fe K band- ~ 39 eV 
at 6.4 keV, a factor of ~ 4 better than the spectral resolution of Suzaku and XMM-Newton. S10 
measured the FWHM of the Fe Ka core, when the data quality was high enough, and found the 
mean FWHM of the Fe Ka core to be 2060 ± 230 km s _1 . When comparing the FWHM of the 
core of Fe K emission to the H/3 optical emission line width, S10 found there to be no universal 
location of the Fe Ka emission relative to the BLR, with possible locations ranging from ~ 0.7-2 
times the radius of the BLR out to the putative torus. 
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We compared the FWHM Fe K core measurements and the corresponding H/3 measurements 
of S10 for those sources common to both S10 and our study. On average, the FWHM measurements 
of the Fe Ka core (6100^320° km s _1 ) of the overlapping sample are comparable to the H/3 FWHM 
(~ 4000 km s _1 ), suggesting that the X-ray line-emitting clouds may be located within the optical 
broad line region. The broader component of Fe Ka emission observed in some objects in our 
sample may ari se from reflection off material from the inner parts of a Compton-thick accretion- 
disk wind (e.g.. lSim et al.ll2008l : lTatum et al.ll2012bl ). The mean a width of the broader component 
was a ~ 350eV, consistent with the findings of Patrick et al. (2012, MNRAS, submitted). 



5.1.2. Spectral variability arguments 



In the context of absorption models, some of the best studied sources hav e had their spectra l 



variability attribu ted to variations in covering- fraction of the absorber (i.e., iTurner et al. 



2008 



Miller et al.ll2008l ). In such a picture, the location of the absorber must be consistent with the 
known time scales for spectral variability. 

Some constraints on gas location have been obtained for sources in our sample. For example, 
our Suzaku observations of NGC 3227 show significant spectral variability on times scales ~ 5 days 
(see the hardness ratio variations in Table [T|), consistent with a virialized cloud location at the 
radius of the BLR. 



In a study of NGC 1365, iRisaliti et al.1 (j2009bl ) compared variations in X-ray hardness, Flux2_5 keV I 
Flux7_io keVi across a 60 ks XMM observation, in an attempt to isolate variations in the column 
of the X-ray absorber. They split the data into three time intervals of ~ 20 ks, based on hardness 
ratio variations. Using these three time intervals in their time-resolved spectroscopy, they found 
that the source occupies two different flux states. The spectral variations were attributed to the 
passage of a cloud having column density Nh ~ 3x 10 23 cm~ 2 crossing the line-of-sight. Assuming 
the cloud was in a Keplerian orbit and cloud density was n = 10 10 cm" 3 , the timescale for the 
event yielded a radial location for the gas r ~ 10 16 cm ( Risaliti et al.ll2009bl ). consistent with the 
radial location and cloud density of the BLR. 

In general, there is irrefutable evidence from X-ray grating spectroscopy that the absorber is 
outflowing (e.g. ITurner et al.ll2008l ). Consequently, the assumption of Keplerian motion is surely 
inaccurate, leading to unreliable estimates of the gas radial location. Independent evidence for 
the location of the X-ray reproce ssor comes from re cent work on X-ray reverberatio n. Estimates 
from sources such as NGC 4051 (jMiller et al.ll2010l ) and Ark 564 (jLegg et al.ll2012l ) suggest that 
reverberation occurs from material tens to hundreds of gravi tational radii fro m the central source. 
However, the reverberating gas may be highly ionized (e.g., iLegg et al.ll2012l ) and may not be the 
same material causing the extreme spectral hardness noted here. 
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5.2. Implications for the intrinsic X-ray luminosity 

With the realization that large columns of gas cover a significant fraction of our line-of-sight, 
we must consider the implications of this Compton-thick absorption to understand the true nuclear 
continuum luminosity. For absorbers with such high column densities as found here, Compton 
scattering has a significant effect and must be taken into account. However, simple estimations of 
the Compton-scattering losses and absorption suppression of the continuum can overestimate the 
effect. The correction factor depends critically on the geometry of the gas, as we need to account 
for light scattered back into the line-of-sight by the cloud complex. 

We can estimate the correction factor between the observed and intrinsic X-ray luminosity in 
the context of the MCRT. For example, the observed 0.5-50 keV luminosity of MCG-03-34-64 (the 
third hardest source in our sample) is L x ~ 1.7 x 10 erg s _1 . After correcting for scattering and 
absorption in the MCRT, the implied intrinsic luminosity is Lj nt (0.5 - 50 keV) ~ 1.7 x 10 44 erg s _1 , 
an order of magnitude higher than that observed. Importantly however, the calculated intrinsi c 



luminosity does not exceed the Eddington luminosity, Lem ~ 10 45 erg s _1 (jMiniutti et al.ll2007bl ). 
suggesting that, in general, partial-covering by a distribution of Compton-thick clouds is consistent 
with the energy budget of the system. 

The correction required to understand the X-ray luminosity does not generally have a signifi- 
cant effect on the bolometric luminosity as the X-ray photons that are hidden from direct view are 
reprocessed into another bandpass and ultimately measured there. 



5.3. Comparison with other hard X-ray selected samples 



Ou r work is not the firs t to determine global properties of AGN using a hard X-ray selected 
sample. IWinter et all (|2009i ) (hereafter W09) examined the X-ray properties of local AGN in the 
BAT 9-month catalog to determine the fraction of "hidden" AGN, which they defined, in part, as 
sources with scattered or leaking light fractions < 3%. Twenty- four percent of AGN in their sample 
were "hidden", most likely by a geometrically-thick torus, and that the type 1 AGN population in 
their sample were relatively unobscured, with low column density gas (log N# < 21.8) along the line- 
of-sight. In our sample, the type 1 AGN show significantly higher fractions of (unabsorbed) leaking 
light (2%-97%) and hence do not meet the W09 criterion for being hidden. Despite not meeting 
the W09 criterion for being hidden, the majority of sources selected in our sample do contain a 
significant amount of Compton-thick, clumpy absorption, a conclusion significantly different to that 
of W09. 



Owing to the differences in our conclusions, we briefly review the most important differences 
in the samples and methodology. The W09 sample includ ed all sources, both type 1 and type 2 
AGN, from the 9-month BAT catalog (jTueller et al.ll2008l ). In contrast, our sample was derived 
from the 58-month BAT catalog, with AGN detection having substantially increased, by about 
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a factor of 4, since the 9-month catalog. Further to this, we cross-correlated the sources in the 
58-month catalog with those in the Suzaku archive and excluded all type 2 AGN, radio-loud AGN, 
and known LINERS. We also note that our AGN type criteria results in ~ 20% more sources (35 in 
W09, 43 in this work) when applied to the 58-month BAT catalog, than if we had applied it to the 
9-month BAT catalog. Even though the 58-month BAT catalog probes deeper flux levels than the 
9-month catalog, our sample is contained in the same observed luminosity parameter space (41.5 < 
log lumi5_50fc e y < 45) as the 9-month study (Figure H|), suggesting that the flux bias in our sample 
selection is negligible. 

There are important methodological differences between the study presented here, and that 
of W09. The spectral analysis of W09 was not restricted to the use of simultaneous data from 
the medium (2 - 10 keV) and hard (> 10 keV) X-ray bands, instead W09 used non-simultaneous 
data from Swift XRT, ASCA SISO, XMM-Newton pn, Chandra ACIS-S, and Suzaku XIS in their 
spectral fitting. Most notably, the spectral fitting performed by W09 was restricted to data below 
10 keV. Only utilizing the bandpass below 10 keV to detect broadband absorption would not be 
sensitive to Compton-thick absorption along the line-of-sight and may result in an underestimation 
of the intrinsic luminosity of the source. W09 compared the absorption-corrected flux of each source 
to the corresponding BAT flux (14 - 195 keV) and found the relationship to be linear, taking this 
as validation of the absorption correction applied to the observed luminosity. However, comparing 
the data below 10 keV to the BAT flux assumes that the time-averaged BAT flux is representative 
of the flux state above 10 keV. Such an assumption is not consistent with th e known variability 



above 10 keV seen in some AGN (e.g. iMiniutti et al.ll2007al : Miller et al.ll2008l ) 



Broadband spectral analysis in the X-ray band offers better constraints on the intrinsic nature 
of AGN, especially for those objects with high column density clouds along the line-of-sight, as the 
model must be self-consistent with data below and above 10 keV. In our analysis, the role of Suzaku 
has been crucial, because it is essential to demonstrate the hard excess with an instrument that 
simultaneously observes the hard and soft X-ray bands: the restriction to use only Suzaku data for 
a BAT-selected sample is a key feature of our study. 



Winter et al.l (12012 . hereafter W12) also selected a sample of sources from the 9-month BAT 
catalog, in order to understand the properties of the AGN warm absorbers. Their sample comprised 
48 type 1 AGN, including both radio-loud and radio-quiet sources up to type 1.5. In this case, 
they used the simultaneous data in the Suzaku archive, where available. For those sources not 
observed by Suzaku, XMM-Newton spectra and time-averaged Swift BAT spectra were used for 
the broadband fits. Their model construction consisted of a powerlaw representing the continuum 
source, a blackbody model parameterizing emission below 2 keV, a Gaussian component modeling 
the Fe Ka emission, and a distant reflector. The type 1 AGN were unobscured, consistent with 
W09, and that the line equivalent width of the Fe Ka emission averaged ~ 100 eV (consistent with 
W09 and our result). The relative strength of the distant reflector spanned 0-5 times that expected 
from a standard thin disk of neutral material subtending 2n sr at the continuum source. Comparing 
the overlapping objects of W12 and our study (excluding all sources with unconstrained relative 
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reflection strength in W12), the average relative strength of the reflector in W12 is 2.10l ' 44 times 
that expected from a standard thin disk of neutral material subtending 2ir sr at the continuum 
source. This result is similar to our observational finding in that the sample is harder than expected, 
especially for a sample of type 1 AGN. The wider span of the relative reflection strength seen in 
the ensemble result of W12 is not only, presumably, a consequence of including both radio-loud and 
radio-quiet AGN in the sample, but also, in part, a consequence of using non-simultaneous data. 
Inclusion of the of radio-loud in a comparison study is problematic, as the nuclear properties and 
the radio jet cannot be separated in the X-ray data. Our study is robust in that we identify general 
properties of a single subset of AGN over a broad energy range using simultaneous data. 



5.4. Comparison to the cosmic X-ray background 



The 1 Ms Chandr a Deep Field South (jGiacconi et all 120021 ). 2 Ms Chandra D eep Field North 
([Alexander et al.ll2003l ) and 800 ks XMM-Newton Lockman Hole (jHasingerl 12004 ) have confirmed 
the prediction of population synthesis models that the cos mic X-ray background originates from 
the integrated emission of extra-galactic point-like sources (|Gillill2004l ). Obscured AGN contribute 
to th e cosmic X-ray backg r ound emission with most of the energ y density being produced above 10 
keV (jComastri et alJl2005l ; iGilli et al.ll2007l ; iTreister et al.ll2009l ). Our observational result suggests 
that a significant fraction of the type 1 AGN population is obscured with at least 50% partial 
covering by Compton-thick material; therefore, it is interesting to compare a model representative 
of our mean spectrum to the spectral shape of the observed cosmic X-ray background. 

As stated earlier, the hardness ratio measured for 1H 0419-577 was consistent with the weighted 
mean of the sample and may represent typical absorber parameters. We assumed 1H 0419-577 was 
repres entative of the sample and used the model construction and parameters from iTurner et al.l 
( 200 9]). We compare d the renormalized cosmic X-ray background data based on Figure 1 of 



Comastri et al. (j2005l ) to the spectral shape of 1H 0419-577 (see Figure [5]). We found the spectral 
shape of 1H 0419-577 to be consistent with the X-ray background above 2 keV with both having 
a characteristic peak ~ 20 - 30 keV. Below 2 keV, the curve for 1H 0419-577 shows a marked rise 
in flux as compared to the X-ray background. These differences are to be expected consid ering 
that emission below 2 keV in type 1 AGN varies with covering fraction ([Miller et al.l 12008) and 
contributions to the X-ray background must include emission from both type 1 and type 2 AGN. 

X-ray population synthesis models typically have source redshifts of about unity. While it 
appears that adding contributions from objects like 1H 0419-577 and shifting them to z = 1, would 
improve the agreement with the X-ray background, any more sophisticated analysis would require 
a proper consideration of the entire local AGN population and its cosmological evolution, which is 
beyond the scope of this paper. 
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6. Conclusion 

We conducted an exploratory study of the hard excess phenomenon in the local type 1 AGN 
population using sources from the Swift BAT 58-month catalog cross-correlated with the holdings 
of the Suzaku public archive. We extracted the hardness ratio and the equivalent width of the total 
Fe Ka emission for each observation and found that a simple disk reflection model does not provide 
an adequate representation of the sample results. Further, the presence of sharp edges in the three 
hardest sources demonstrate that the extreme hardness cannot be dominated by blurred reflection 
components. 

However, a model comprising partial-covering by an ensemble of Compton-thick clouds is 
consistent with the observational result. In the context of such a model, ~ 80% of the sample had a 
hardness ratio consistent with > 50% covering of the continuum by low-ionization, Compton-thick 
gas. Simple considerations suggested that the absorber lies at or within the optical BLR. Future 
work will include the results of our detailed analysis of this sample and a comparison of the results 
of this work with model predictions from a ionized, Compton-thick, partial-covering absorber, that 
includes Compton-scattering. More extensive spectral analysis using a parameterization of the 
Monte-Carlo radiative transfer model can further restrict the gas location as well as offer a better 
understanding of the geometry of the ensemble of Compton-thick clouds. 

7. Acknowledgements 

MMT thanks Adam Patrick and Jason Gofford for aiding in the data reduction. MMT would 
like to acknowledge NASA grant NNX08AJ41G. TJT would like to acknowledge NASA grant 
NNX11AJ57G. LM acknowledges support from STFC grant ST/H002456/1. 

REFERENCES 

Alexander, D. M., Bauer, F. E., Brandt, W. N., Schneider, D. P., Hornschemeier, A. E., Vignali, 
C, Barger, A. J., Broos, P. S., Cowie, L. L., Garmire, G. P., Townsley, L. K., Bautz, M. W., 
Chartas, G., k Sargent, W. L. W. 2003, AJ, 126, 539 

Anders, E. & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 

Baumgartner, W. H., Tueller, J., Markwardt, C, & Skinner, G. 2010, in Bulletin of the American 
Astronomical Society, Vol. 42, AAS/High Energy Astrophysics Division #11, 675 

Blustin, A. J., Kriss, G. A., Holczer, T., Behar, E., Kaastra, J. S., Page, M. J., Kaspi, S., 
Branduardi-Raymont, G., & Steenbrugge, K. C. 2007, A&A, 466, 107 

Blustin, A. J., Page, M. J., Fuerst, S. V., Branduardi-Raymont, G., & Ashton, C. E. 2005, A&A, 
431, 111 



-19- 



Boldt, E. 1987, in IAU Symposium, Vol. 124, Observational Cosmology, ed. A. Hewitt, G. Burbidge, 
k L. Z. Fang, 611-615 

Brightman, M. k Nandra, K. 2012, arXivl202.1291 

Comastri, A., Gilli, R., k Hasinger, G. 2005, Experimental Astronomy, 20, 41 
Comastri, A., Setti, G., Zamorani, G., k Hasinger, G. 1995, A&A, 296, 1 
Crenshaw, D. M. k Kraemer, S. B. 1999, ApJ, 521, 572 
Crenshaw, D. M., Kraemer, S. B., k George, I. M. 2003, ARA&A, 41, 117 
Dickey, J. M. k Lockman, F. J. 1990, ARA&A, 28, 215 

Ebisawa, K., Ogawa, M., Aoki, T., Dotani, T., Takizawa, M., Tanaka, Y., Yoshida, K., Miyamoto, 
S., Iga, S., Hayashida, K., Kitamoto, S., k Terada, K. 1994, PASJ, 46, 375 

George, I. M. k Fabian, A. C. 1991, MNRAS, 249, 352 

Giacconi, R., Zirm, A., Wang, J., Rosati, P., Nonino, M., Tozzi, P., Gilli, R., Mainieri, V., Hasinger, 
G., Kewley, L., Bergeron, J., Borgani, S., Gilmozzi, R., Grogin, N., Koekemoer, A., Schreier, 

E. , Zheng, W., k Norman, C. 2002, ApJS, 139, 369 

Gilli, R. 2004, Advances in Space Research, 34, 2470 
Gilli, R., Comastri, A., k Hasinger, G. 2007, A&A, 463, 79 

Gilli, R., Vignali, C, Mignoli, M., Iwasawa, K., Comastri, A., k Zamorani, G. 2010, A&A, 519, 
A92 

Haardt, F. k Maraschi, L. 1993, ApJ, 413, 507 

Hasinger, G. 2004, Nuclear Physics B Proceedings Supplements, 132, 86 

Heckman, T. M., Ptak, A., Hornschemeier, A., k Kauffmann, G. 2005, ApJ, 634, 161 

Kallman, T. k Bautista, M. 2001, ApJS, 133, 221 

Kaspi, S., Brandt, W. N., George, I. M., Netzer, H., Crenshaw, D. M., Gabel, J. R., Hamann, 

F. W., Kaiser, M. E., Koratkar, A., Kraemer, S. B., Kriss, G. A., Mathur, S., Mushotzky, 
R. F., Nandra, K., Peterson, B. M., Shields, J. C, Turner, T. J., k Zheng, W. 2002, ApJ, 
574, 643 

Krolik, J. H. k Kallman, T. R. 1987, ApJ, 320, L5 

LaMassa, S. M., Heckman, T. M., Ptak, A., Martins, L., Wild, V., Sonnentrucker, P., k Horn- 
schemeier, A. 2011, ApJ, 729, 52 



-20- 



Lee, J. C, Ogle, P. M., Canizares, C. R., Marshall, H. L., Schulz, N. S., Morales, R., Fabian, A. C, 
& Iwasawa, K. 2001, ApJ, 554, L13 

Legg, E., Miller, L., Turner, T. J., Giustini, M., Reeves, J. N., & Kraemer, S. B. 2012, 
arXiv: 1210.0469 

Maeda, Y., Someya, K., Ishida, M., the XRT team, Hayashida, K., Mori, H., & the XIS team. 
2008, JX-ISAS-SUZAKU-MEMO-2008-06 

Magdziarz, P. & Zdziarski, A. A. 1995, MNRAS, 273, 837 

McKernan, B. k Yaqoob, T. 1998, ApJ, 501, L29 

McKernan, B., Yaqoob, T., & Reynolds, C. S. 2007, MNRAS, 379, 1359 
Miller, L., Turner, T. J., & Reeves, J. N. 2008, A&A, 483, 437 
-. 2009, MNRAS, 399, L69 

Miller, L., Turner, T. J., Reeves, J. N., George, I. M., Kraemer, S. B., & Wingert, B. 2007, A&A, 
463, 131 

Miller, L., Turner, T. J., Reeves, J. N., Lobban, A., Kraemer, S. B., & Crenshaw, D. M. 2010, 
MNRAS, 403, 196 

Miniutti, G., Fabian, A. C, Anabuki, N., Crummy, J., Fukazawa, Y., Gallo, L., Haba, Y., 
Hayashida, K., Holt, S., Kunieda, H., Larsson, J., Markowitz, A., Matsumoto, C, Ohno, 
M., Reeves, J. N., Takahashi, T., Tanaka, Y., Terashima, Y., Torii, K., Ueda, Y., Ushio, M., 
Watanabe, S., Yamauchi, M., & Yaqoob, T. 2007a, PASJ, 59, 315 

Miniutti, G., Ponti, G., Dadina, M., Cappi, M., & Malaguti, G. 2007b, MNRAS, 375, 227 

Murphy, K. D. & Yaqoob, T. 2009, MNRAS, 397, 1549 

Nandra, K. & George, I. M. 1994, MNRAS, 267, 974 

Nandra, K., O'Neill, P. M., George, I. M., & Reeves, J. N. 2007, MNRAS, 382, 194 

Netzer, H., Kaspi, S., Behar, E., Brandt, W. N., Chelouche, D., George, I. M., Crenshaw, D. M., 
Gabel, J. R., Hamann, F. W., Kraemer, S. B., Kriss, G. A., Nandra, K., Peterson, B. M., 
Shields, J. C, & Turner, T. J. 2003, ApJ, 599, 933 

Reeves, J. N., Nandra, K., George, I. M., Pounds, K. A., Turner, T. J., & Yaqoob, T. 2004, ApJ, 
602, 648 

Reeves, J. N., O'Brien, P. T., Braito, V., Behar, E., Miller, L., Turner, T. J., Fabian, A. C, Kaspi, 
S., Mushotzky, R., & Ward, M. 2009, ApJ, 701, 493 



- 21 - 



Risaliti, G., Bianchi, S., Matt, G., Baldi, A., Elvis, M., Fabbiano, G., & Zezas, A. 2005, ApJ, 630, 
L129 

Risaliti, G., Braito, V., Laparola, V., Bianchi, S., Elvis, M., Fabbiano, G., Maiolino, R., Matt, G., 
Reeves, J., Salvati, M., & Wang, J. 2009a, ApJ, 705, LI 

Risaliti, G., Miniutti, G., Elvis, M., Fabbiano, G., Salvati, M., Baldi, A., Braito, V., Bianchi, S., 
Matt, G., Reeves, J., Soria, R., & Zezas, A. 2009b, ApJ, 696, 160 

Scott, A. E., Stewart, G. C., Mateos, S., Alexander, D. M., Hutton, S., & Ward, M. J. 2011, 
MNRAS, 417, 992 

Shu, X. W., Yaqoob, T., & Wang, J. X. 2010, ApJS, 187, 581 

Sim, S. A., Long, K. S., Miller, L., & Turner, T. J. 2008, MNRAS, 388, 611 

Sim, S. A., Miller, L., Long, K. S., Turner, T. J., & Reeves, J. N. 2010, MNRAS, 404, 1369 

Tanaka, Y., Boiler, T., Gallo, L., Keil, R., & Ueda, Y. 2004, PAS J, 56, L9 

Tatum, M. M., Turner, T. J., Miller, L., & Reeves, J. N. 2012a, in prep. 

Tatum, M. M., Turner, T. J., Sim, S. A., Miller, L., Reeves, J. N., Patrick, A. R., & Long, K. S. 
2012b, ApJ, 752, 94 

Tombesi, F., Cappi, M., Reeves, J. N., Palumbo, G. G. C., Braito, V., & Dadina, M. 2011, ArXiv 
e-prints 

Tombesi, F., Cappi, M., Reeves, J. N., Palumbo, G. G. C., Yaqoob, T., Braito, V., & Dadina, M. 
2010, A&A, 521, A57 

Treister, E., Urry, C. M., & Virani, S. 2009, ApJ, 696, 110 

Tueller, J., Mushotzky, R. F., Barthelmy, S., Cannizzo, J. K., Gehrels, N., Markwardt, C. B., 
Skinner, G. K., & Winter, L. M. 2008, ApJ, 681, 113 

Turner, T. J. & Miller, L. 2009, A&A Rev., 17, 47 

Turner, T. J., Miller, L., Kraemer, S. B., & Reeves, J. N. 2011, ApJ, 733, 48 

Turner, T. J., Miller, L., Kraemer, S. B., Reeves, J. N., k Pounds, K. A. 2009, ApJ, 698, 99 

Turner, T. J., Miller, L., Reeves, J. N., k Kraemer, S. B. 2007, A&A, 475, 121 

Turner, T. J., Reeves, J. N., Kraemer, S. B., & Miller, L. 2008, A&A, 483, 161 

Veron-Cetty, M.-P. & Veron, P. 2006, A&A, 455, 773 

Walton, D. J., Reis, R. C., & Fabian, A. C. 2010, MNRAS, 408, 601 



- 22 - 



Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914 

Winter, L. M., Mushotzky, R. F., Reynolds, C. S., k Tueller, J. 2009, ApJ, 690, 1322 
Winter, L. M., Veilleux, S., McKernan, B., & Kallman, T. R. 2012, ApJ, 745, 107 
Yaqoob, T., Murphy, K. D., Miller, L., & Turner, T. J. 2010, MNRAS, 401, 411 



This preprint was prepared with the AAS IATjrjX macros v5.2. 



- 23 - 



o- 



I 

X 



o 

I 

X 



+ 
+ 



X 



| fl it 



— *t 



A 



10 20 



50 100 200 



500 



BAT Flux (10 12 ergs s 1 cm- 2 ) 



Fig. 1. — : The hardness ratio, Fluxi5_50fcev/Flux2_iofc e y, plotted against the Swift BAT flux. 
Superimposed is the weighted hardness ratio mean (solid horizontal line). Overlaid are model 
predictions for simple disk reflection for R=l (dashed line) and R— > oo (dash-dot) and for partial 
covering by neutral, Compton-thick, N#=2 x 10 24 cm' 2 , clouds with 98% (circle-dot), 90% (circle- 
cross), 70% (square), and 50% (triangle) covering fractions, not accounting for Compton-scattering 
losses. 
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Fig. 2. — The hardness ratio plotted against the total Fe Ka equivalent width. (Left) Model 
lines are overlaid for disk reflection from R=l to pure reflection (thick dash light blue); for a 
Compton-thick torodial reprocessor at an inclination of 60° (small thick solid gray) and 70° (large 
thick solid gray); and for a thin absorbing shell (thick dash-dot orange). MCRT models comprise 
the loops (solid red, dash green, dash-dot blue, dot cyan, dash-dot-dot-dot magenta). (Right) 
MCRT model predictions in more detail: each closed loop corresponds to a particular distribution 
of clouds. Moving counterclockwise around any one loop corresponds to increasing column density 
of an individual spherical cloud. See text for details 
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Fig. 3. — Data/model ratio residuals from Suzaku observations of NGC 1194 (top left), NGC 1365 
(2010 July 15, top right) and MCG-3-34-64 (below left). For each observation, the 5-7 keV data 
were fit with a powerlaw and a Gaussian component at 6.4 keV. The model was extrapolated to 9 
keV, and the 5-9 keV data were divided by the folded model. Each spectrum shows a sharp edge 
at ~ 7 keV. 
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Fig. 4. — The 15-50 keV luminosity distribution for those sources in our sample that are contained 
in the 9-month (solid black) and 58-month (dash red) BAT catalog. 
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Fig. 5. — The best -fitting model for 1H 0419-577 (solid line), based on the model parameters from 
Turner et al.l (|2009i ). compared with t he ren ormalized cosmic X-ray background data (dashed line) 
based on Figure 1 of lComastri et al.l (| 20051 ). 
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Table 1. X-ray Observation Source List 



Object 


ObsID 


RA (h 


m s) 1 


t-. / r\ I II \ 1 

Dec (° ' ") 


1 

z 1 


N/J (Gal) 2 


m 1 

Type 1 


Hardness Ratio 3 


Total EW (eV) 4 


1H 0419-577 


702041010 


04 


26 


00.7 


-57 12 01 


0.104 


0.020 


1.5 


1.45+0.06 


42+ 1 

— 1 


1H 0419-577 


704064010 


04 


26 


00.7 


-57 12 01 


0.104 


0.020 


1.5 


1.53+0.07 


20+1 


2MASXJ09043699+5536025 


704027010 


09 


04 


36.9 


+55 36 03 


0.037 


0.028 


1 


3.24+2.30 


97 + r 

" —5 


Ark 120 


702014010 5 


05 


16 


11.4 


-00 08 59 


0.033 


0.126 


1 


1.10+0.05 


197+} 


ESO 362-18 


703014010 


05 


19 


35.8 


-32 39 28 


0.012 


0.018 


1.5 


1.44+0.55 


175+i 


ESO 548-G081 


704026010 5 


03 


42 


03.7 


-21 14 39 


0.015 


0.032 


l 6 


1.66+0.14 


94+? 
— l 


Fairall 9 


702043010 5 


01 


23 


45.8 


-58 48 21 


0.047 


0.032 


1.2 


1.31+0.04 


186+1 


Fairall 9 


705063010 


01 


23 


45.8 


-58 48 21 


0.047 


0.032 


1.2 


1.58+0.04 


78+2 


IC4329A 


702113010 5 


13 


49 


19.3 


-30 18 34 


0.016 


0.044 


1.2 


1.57+0.04 


115+1! 

— 15 


IC4329A 


702113030 


13 


49 


19.3 


-30 18 34 


0.016 


0.044 


1.2 


1.66+0.04 


58+"' 

° — 10 


IC4329A 


702113040 5 


13 


49 


19.3 


-30 18 34 


0.016 


0.044 


1.2 


1.49+0.04 


126+1 


IC4329A 


702113050 


13 


49 


19.3 


-30 18 34 


0.016 


0.044 


1.2 


1.90+0.05 


98+" 

^ — 19 


IC4329A 


702113020 


13 


49 


19.3 


-30 18 34 


0.016 


0.044 


1.2 


1.43+0.03 


50+;? 

— 11 


ICR J16185-5928 


702123010 


16 


18 


36.4 


-59 27 17 


0.035 


0.247 


1 NL 6 


1.67+0.15 


75+?," 

' — 22 


IGR J19378-0617 


705055010 


19 


37 


33.0 


-06 13 05 


0.010 


0.148 


1.5 


0.68+0.07 


47+;? 

1 -15 


MCG+04-22-42 


704028010 


09 


23 


43.0 


+22 54 33 


0.032 


0.034 


1.2 


1.27+0.14 


95+2 


MCG+08-11-11 


702112010 


05 


54 


53.6 


+46 26 22 


0.021 


0.209 


1.5 


1.54+0.03 


100+?,, 

— 10 


MCG-02-14-09 


703060010 5 


05 


16 


21.2 


-10 33 41 


0.029 


0.093 


1 


1.44+0.25 


200+2 


MCG-02-58-22 


704032010 


23 


04 


43.5 


-08 41 09 


0.047 


0.036 


1.5 


1.82+0.03 


47+5 

' —6 


MCG-03-34-64 


702112010 


13 


22 


24.4 


-16 34 42 


0.017 


0.058 


1 h 


7.53+0.56 


372+!!5 


MCG-06-30-15 


100004010 


13 


35 


53.7 


-34 17 44 


0.008 


0.041 


1.2 


1.24+0.05 


98+J 


MCG-06-30-15 


700007010 


13 


35 


53.7 


-34 17 44 


0.008 


0.041 


1.2 


1.33+0.03 


80+" 


MCG-06-30-15 


700007020 


13 


35 


53.7 


-34 17 44 


0.008 


0.041 


1.2 


1.65+0.04 


63+g 1 


MCG-06-30-15 


700007030 


13 


35 


53.7 


-34 17 44 


0.008 


0.041 


1.2 


1.59+0.03 


49+1 


MR 2251-178 


704055010 


22 


54 


05.8 


-17 34 55 


0.064 


0.027 


1.5 


1.34+0.04 


38_ 1 3 


Mrk 1018 


704044010 


02 


06 


16.0 


-00 17 29 


0.042 


0.026 


1.9 


1.73+0.14 


76+2 


Mrk 110 


702124010 


09 


25 


12.9 


+52 17 11 


0.035 


0.015 


1 


1.36+0.06 


64+; 3 

" — 18 


Mrk 279 


704031010 


13 


53 


03.4 


+69 18 30 


0.031 


0.018 


1 


3.21+0.21 


234+3 


Mrk 335 


701031010 5 


00 


06 


19.5 


+20 12 10 


0.026 


0.040 


1 


1.05+0.09 


181+!* 

° — 12 


Mrk 352 


704025010 


00 


59 


53.3 


+31 49 37 


0.015 


0.055 


1 


1.42+0.28 


89+3? 


Mrk 359 


701082010 


01 


27 


32.5 


+19 10 44 


0.017 


0.048 


1 NL 


1.38+0.44 


1351^ 


Mrk 509 


701093010 


20 


44 


09.7 


-10 43 25 


0.034 


0.041 


1.5 


1.41+0.06 




Mrk 509 


701093020 


20 


44 


09.7 


-10 43 25 


0.034 


0.041 


1.5 


1.56+0.06 


sell? 


Mrk 509 


701093030 


20 


44 


09.7 


-10 43 25 


0.034 


0.041 


1.5 


1.51+0.08 




Mrk 509 


701093040 


20 


44 


09.7 


-10 43 25 


0.034 


0.041 


1.5 


1.50+0.06 




Mrk 766 


701035010 


12 


18 


26.5 


+29 48 46 


0.013 


0.017 


1 NL 


1.27+0.11 


DJ -33 


Mrk 766 


701035020 


12 


18 


26.5 


+29 48 46 


0.013 


0.017 


1 NL 


1.18+0.10 




Mrk 79 


702044010 


07 


42 


32.8 


+49 48 35 


0.022 


0.053 


1.2 


1.78+0.16 


15911^ 


Mrk 841 


701084010 


15 


04 


01.2 


+10 26 16 


0.036 


0.023 


1.5 


1.22+0.09 


83+1 


Mrk 841 


701084020 


15 


04 


01.2 


+10 26 16 


0.036 


0.023 


1.5 


2.03+0.12 


58l?° 


NGC 1194 


704043010 


03 


03 


49.1 


-01 06 13 


0.014 


0.071 


1.9 


31.17+2.92 


789+121 
' sa -140 


NGC 1365 


702047010 


03 


33 


36.4 


-36 08 25 


0.005 


0.014 


1.8 


3.07+0.08 


33ll« 


NGC 1365 


705031010 


03 


33 


36.4 


-36 08 25 


0.005 


0.014 


1.8 


11.12+0.33 




NGC 1365 


705031020 


03 


33 


36.4 


-36 08 25 


0.005 


0.014 


1.8 


12.72+0.22 


52311 
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Table 1 — Continued 



Object 


ObsID 


RA (h m s) 1 


Dec (° ' "f 


z 1 




(Gal) 2 


Type 1 


Hardness Ratio 3 


Total EW (eV) 4 


INCjL' ozz ( 


I U3U22010 


10 


23 30.6 


+ ly 51 54 


0.004 





.020 


1 . 5 


l.y5±0.43 


,,,+ 14 

113_ 17 


INLjU oil 1 


no on on 
I 03022020 


10 


23 30.6 


i i n ci c a 
+ 1» 51 54 


0.004 


0. 


020 


1.5 


n n7 I n 1 o 

3. y (±0.12 


r, fi C + 20 

z»O_20 


tvt goo7 

INLjL/ Qllt 


*7nQnoonQn 
(03022030 


10 


23 30.6 


+ 1» 51 54 


0.004 


0. 


.020 


1.5 


2.0 ( ±0.08 


163_ 19 


INLiO ozz/ 


1 03022040 


10 


23 30.6 


i 1 n ci c a 

+iy 51 54 


0.004 


0. 


020 


1.5 


O.2o±0.22 


497_ 34 


AT (~1 (~\ Q007 


/03022050 


10 


23 30.6 


I in ci er a 

-|-iy 51 54 


0.004 


0. 


.020 


1.5 


o noxfi m n 
2.y2±0.10 


162l^ 3 


INLiL/ OZZI 


vnqnoon^n 
( 030220b0 


10 


23 30.6 


i 1 n c 1 c a 

+iy 51 54 


0.004 


0. 


020 


1.5 


2.6/±0.11 


lb3_ 17 


INLjO o51b 


100031010 


11 


06 47.5 


1 TO OA nT 

-f-72 o4 U/ 


0.009 


0. 


.031 


1.5 


2.y/±0.0o 


215_ 12 


INCjU oolb 


(04002010 


11 


06 47.5 


1 TO Q A nT 

+ J Z o4 U/ 


0.009 


0. 


031 


1.5 


2.58±0.0 ( 


no i ~r 2t) 

281_ 19 


INLjU o/oo 


*7m no on 1 n 
(01033010 


11 


39 01.7 


Q7 A A in 
— o7 44 iy 


0.010 


0. 


.085 


1.5 


2.02±0.0o 


164_ 10 


INLtL/ o/oo 


/ U4U0OU1U 


11 


39 01.7 


07 /i ^ in 
— ol 44 iy 


n ni n 
U.U1U 


n 
U. 


Uo5 


1.0 


l.yb±U.Uo 


lU2_g 


mpp /inKi 
INLiL/ 4U51 


(UUUU4U1U 


12 


03 09.6 


+44 ol 5o 


n nno 
U.UUZ 


n 
U. 


Ulo 


1 

1 


O *7KJ-C\ 1 n 

Z. i o±U. iy 


1S1+18 
181_ 13 


INLjO 4U51 


vnqnoQm n 


12 


03 09.6 


+44 ol 5o 


n nno 
U.UUz 


n 
U 


Ulo 


1 

1 


1 A 1 J-Ci HA 

1.41±U.U4 


»2_ 10 


INLjO 4U51 


vnqnoqnon 


12 


03 09.6 


i a a q i c;q 
+44 ol 5o 


n nno 
U.UUz 


n 
U. 


Ulo 


1 

1 


i.oy±u.iu 


m+18 


INLjL; 4151 


/ UlUo4UlU 


12 


10 32.6 


1 QO OA 01 

+oy z4 zi 


n nnQ 
U.UUo 


n 
U. 


non 
UzU 


"1 K 
1.0 


Q TKJ-(\ Oft 
O. / 0±U.UD 


OOQ+19 




(UZU4UU1U 


12 


39 39.4 


nc; on Qn 
— UO zU oy 


n nnn 
U.UUy 


n 
U. 


noQ 
UZo 


1 

1 


l.yo±U.14 


278_ 21 


Mfip re nft 
INLjL/ 55UO 


7m nqnm n5 
/ UluoUUlU 


14 


13 14.9 


HQ 1 O 07 

— Uo 1Z Zi 


n nnft 
U.UUO 


n 
U. 


nog 

Uoo 


1 l 


i.t>y±u.uo 


1,9+7 


mpp nft 


(UlUoUUZU 


14 


13 14.9 


HQ 1 O 07 

— Uo 1Z z / 


n nnft 
U.UUO 


n 
U. 


nQC 
Uoo 


1 l 


1. i U±U.Uo 


i on+7 
1ZU ij. 


NGC 5506 


701030030 


14 


13 14.9 


-03 12 27 


0.006 


0. 


.038 


1 i 


1.69±0.03 


75l?o 


NGC 5548 


702042010 


14 


17 59.5 


+25 08 12 


0.017 


0. 


017 


1.5 


3.52±0.34 


193+^ 


NGC 5548 


702042020 


14 


17 59.5 


+25 08 12 


0.017 


0. 


.017 


1.5 


2.29±0.16 


151±1 


NGC 5548 


702042040 


14 


17 59.5 


+25 08 12 


0.017 


0. 


017 


1.5 


1.57±0.10 


o 9 +19 


NGC 5548 


702042050 


14 


17 59.5 


+25 08 12 


0.017 


0. 


.017 


1.5 


2.24±0.14 


OQ+14 
yo_ 12 


NGC 5548 


702042060 


14 


17 59.5 


+25 08 12 


0.017 


0. 


.017 


1.5 


1.70±0.08 




NGC 5548 


702042070 


14 


17 59.5 


+25 08 12 


0.017 


0. 


017 


1.5 


2.03±0.11 




NGC 5548 


702042080 


14 


17 59.5 


+25 08 12 


0.017 


0. 


.017 


1.5 


3.08±0.21 


15611 


NGC 6860 


703015010 s 


20 


08 46.9 


-61 06 01 


0.015 


0. 


.042 


1.5 


1.42±0.08 


106 + } 


NGC 7213 


701029010 


22 


09 16.3 


-47 10 00 


0.006 


0. 


.020 


1.5 6 


1.44±0.06 


891" 


NGC 7314 


702015010 


22 


35 46.2 


-26 03 02 


0.005 


0. 


.015 


1 h 


1.56±0.15 


138lg 


NGC 7469 


703028010 


23 


03 15.6 


+08 52 26 


0.016 


0. 


.049 


1.5 


1.83±0.06 


149li 3 6 


NGC 985 


704042010 


02 


34 37.8 


-08 47 15 


0.043 


0. 


.029 


l 6 


1.65±0.17 




SWIFT J2127.4+5654 


702122010 5 


21 


27 44.9 


+56 56 40 


0.014 


0. 


.787 


1 


0.90±0.04 


159+1° 


UGC 06728 


704029010 


11 


45 16.0 


+79 40 53 


0.007 


0. 


.045 


1.2 


1.28±0.12 


OQ + 31 

»3_ 29 



^Veron-Cettv fc Veronl fcOOcj l 

2 The Galacti c column density in units of 10 22 cm -2 obtained from the weighted average Njj in the Dickey &; Lockman HI in the 
Galaxy survey l lDickev fc Lockmaniriggd ) 

3 Flux 15 _5o fc e v/F mx 2-io feeV with errors calculated from the net count rate errors of the XIS and PIN data 

4 Total equivalent width of the Fe Ka emission line with the errors to 90% confidence 

5 Required two Gaussian components to parameterize the Fe Ka emission. The EW for these observations were the sum of the narrow 
and broad base of the Fe K emission 

6 Sp ectral type as defined in the Swift BAT 58-month source catalog. Used when no spectral type was available in lVeron-Cettv &: Veronl 
j2006h . 



